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^ , Abstract 

>: 

^ ■ QCD Laplace transform sum rules, involving the axial-vector current divergences, 

are used in order to determine the strange quark mass. The two-point function is 
known in QCD up to four loops in perturbation theory, and up to dimension-six in 
the non-perturbative sector. The hadronic spectral function is reconstructed using 
threshold normalization from chiral symmetry, together with experimental data for 
the two radial excitations of the kaon. The result for the running strange quark mass, 
in the WS scheme at a scale of 1 GeV^ is: rhsilGeV'^) = 155 ± 25MeV. 



In spite of many attempts |T[]-[§] to improve the accuracy of QCD sum rule determinations 
of the strange quark mass, present uncertainties still remain uncomfortably large. This is a 
serious limiting factor affecting areas such as kaon-physics and CP violation, which depend 
strongly on m^. One of the determinations believed to be among the most accurate was 
the one that combined the current algebra ratio ^ 

TfX 

12.6 ±0.5, (1) 



rriu + md 

together with a QCD sum rule determination of (m^ + rri^). The latter was discussed some 
years ago in the framework of QCD Finite Energy Sum Rules (FESR) in the pseudoscalar 
channel, at the two-loop level in perturbative QCD, and including non-perturbative con- 
densates up to dimension-six, with the result 

(m„ + md) (1 GeV) = 15.5 ±2.0 MeV. (2) 

This result, together with Eq.(l), implies 

rus (1 GeV) = 195 ± 28 MeV. (3) 

A recent re-analysis ^ of the same QCD sum rules, but including the next (3-loop) order 
in perturbation theory, and a somewhat different hadronic spectral function, obtains 

{rriu + rud) (1 GeV) = 12.0 ± 2.5 MeV , (4) 

which, using Eq.(l), leads to 

m, (1 GeV) = 151 ± 32 MeV . (5) 

The same raw data for resonance masses and widths, plus the same threshold normalization 
from chiral perturbation theory, has been used in both analyses [Q and §]. The difference 
in the results cannot be accounted for by the inclusion (or not) of the 3-loop perturbative 
QCD contribution, which amounts to a reduction of the 2-loop result, Eq.(2), of only a 
few percent. Instead, this difference stems mainly from the different reconstructions of the 
spectral function from the resonance data; with a more elaborate functional form being 



adopted in [Q. This reveals a type of systematic uncertainty of the QCD sum rule method, 
that will cease to be an uncertainty only after the pseudoscalar hadronic spectral function 
is measured directly and accurately (e.g. from tau-lepton decays). In the absence of such 
direct data, it is not possible to exclude the presence of a background (constructive or de- 
structive) beyond that implicit in the chiral normalization. Also, this overall normalization 



could assume a non-standard value as advocated in |Ty]. Hence, at present, both results 
for rris, Eqs. (3) and (5), are equally acceptable, and taken together provide a measure of 
underlying systematic uncertainties. 

A second alternative is based on QCD sum rules involving the correlator of the strange 
vector current divergences. In this case, the availability of experimental data on K — tt 
phase shifts [^ allows, in principle, for a reconstruction of the hadronic spectral function 
in this channel, from threshold up to s ~ 7 GeV . In both [^] and [^, the functional form 
chosen for this reconstruction consisted of a superposition of two Breit-Wigner resonances, 
corresponding to the Kq (1430) and the Kq (1950) |12[, normalized at threshold according 



to conventional chiral-symmetry. It was argued in [Q and that the non-resonant back- 
ground implicit in this threshold normalization was important to achieve a good fit to the 
K ~ 7T phase shifts. In |^ and |^, the relevant correlator was calculated in perturbative 
QCD at the 3-loop level, with mass corrections up to the quartic order, and including 
non-perturbative quark and gluon vacuum condensates up to dimension four (the d = 6 
condensates are numerically unimportant and can be safely neglected ||5|). Using Laplace 
transform sum rules, the results for the strange quark mass thus obtained were 

171±15MeV ([4]) 
m, (1 GeV) = (6) 

178±18MeV (§) 

The errors reflect uncertainties in the experimental data, in the value of Aqqj) (Aqqj) ~ 
200 - 500 MeV), in the continuum threshold sq (sq — 6 - 7 GeV^), and in the values of the 
vacuum condensates. 

As mentioned earlier, an example of a potential systematic error affecting these results 
would be the presence of a background, beyond the one implicit in the chiral-symmetry 



normalization of the hadronic spectral function at threshold. Obviously, this is not included 
in (6). A reanalysis of this QCD sum rule determination of rus has uncovered this 
uncertainty. In fact, it is claimed in ^ that by using the Omnes representation to relate 
the spectral function to the K — n phase shifts, it is necessary to include a background 
interfering destructively with the resonances. As a result, the hadronic spectral function is 
considerably smaller than that used in 0-0. This in turn implies smaller values of m^, 
viz. 

m, (1 GeV) = 140 ± 20MeV . (7) 

Still another source of systematic uncertainty, this time of a theoretical nature, has been 
unveiled in [^. This has to do with the following two possibilities: (a) expand the QCD 
correlator in inverse logarithms of the momentum transfer Q^, or (b) expand only in terms 
of powers of the strong coupling a^. Similarly, after Laplace transforming the correlator 
one faces the same problem, except that the momentum transfer is being replaced by the 
Laplace variable M^. It has been argued in P] that it makes more sense to make full use of 
the perturbative expansions of the quark mass and coupling (known to 4-loop order), and 
hence not to expand them. Numerically, it turns out that the non-expanded expression 
is far more stable than the truncated one, when moving from one order in perturbation 
theory to the next. This fact lends strong support to the non-expanded alternative. In 
addition, as shown in P], logarithmic truncation can lead to sizable overestimates of ra- 
diative corrections. This in turn implies an underestimate of the quark mass. In fact, 
after using untruncated expressions, together with the same hadronic spectral function 
parametrization as in [@]-[H, the authors of P] find 

m, (1 GeV) = 203 ± 20MeV , (8) 

to be compared with the results (6) obtained from truncated expressions. Until the ques- 
tions of truncation, and of the correct form of the hadronic spectral function become 
satisfactorily settled, one should take the value (8) together with (6), and include (7) as 
well. This would give nis (1 GeV) = 170 ± 50MeV, a rather inaccurate, albeit realistic 



result. 



In this note we discuss a direct determination of rris using QCD Laplace sum rules in the 
pseudoscalar channel, i.e. involving the correlator 



Mq' 



I (fx e'^" < Q\T{d^ A^{x) d" 4(0)) |0 > 



(9) 



where A^{x) =: s(x)7^75-u(x) :, and d^ A^{x) = m^ : s(x)i75 u{x) : . The QCD expression 
of this two-point function is known |Q,|^, at the four-loop level in perturbative QCD, 
and up to dimension six in the non-perturbative sector. Also, the old problem of mass 
singularities has been satisfactorily solved in m,[§]. As a result of this, quark mass correc- 
tions are also known up to quartic order. Notice that the QCD result for the correlator (9) 
is trivially obtained from that involving the vector current divergences; hence, they both 
look quite similar. The QCD expression of the Laplace transform of Eq.(9), i.e. 



4iM' 



4iQ' 



1 



e-'/^'^ - Im V'5(s) ds 
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is given by 
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and where 7^ is Euler's constant, ({n) is Riemann's zeta function, oi = 2795.0778, all 
numerical coefficients refer to three flavours and three colours, and we have neglected the 
up-quark mass everywhere. Given the uncertainties of the method, plus the size of sys- 
tematic errors, it is not justifled to keep m„ different from zero. The four- loop expressions 
for the strong running coupling and quark mass are given by 
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where L = IoqIM^/Aqq,^), LL = logL, and 



ft._!Hi£?-™c(3). 



4608 



32 



(17) 



with 74 = 88.5258 (see [0). In line with the discussion after Eq. (7), and following 0, we 
shall not expand the above QCD expressions in inverse powers of L, but rather substitute 
the numerical values of ^^(M^) and rhslM"^) as determined from Eqs.(15)-(16) for a given 
value of Aqcd- The dimension-six non-perturbative term has been omitted as it is of no 
numerical importance. 



The hadronic spectral function associated with the correlator (9) is very different from that 

of the vector divergences. There is, at present, preliminary information from tau-decays 

1^ in a kinematical range restricted by the tau-mass. We reconstruct the spectral func- 



tion, including in addition to the kaon-pole its radial excitations K(1460) and K(1830), 
normalized at threshold according to conventional chiral symmetry. In addition, we incor- 
porate the resonant sub-channel i^*(892) — vr, which is of numerical importance given the 
narrow width of the i^*(892) (the sub-channel p(770) — K is numerically negligible). This 
chiral symmetry normalization is of the form 
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The pion mass has been neglected above, in line with the approximation ttIu = made in 
the QCD sector, and in our normalization /,r — 93MeV". The complete hadronic spectral 
function is then 



TV 71 [i -\- A) 

+ -lm^5{s)\QCD0{s-so) , (21) 

n 

where fx — ^■'^^fn, Im 'ip5{s)\QCD is the perturbative QCD spectral function modelhng the 
continuum which starts at some threshold Sq, BWi^2{s) are Breit-Wigner forms for the two 
kaon radial excitations, normalized to unity at threshold, and A controls the relative im- 
portance of the second radial excitation. The choice A ~ 1 results in a reasonable (smaller) 
weight of the K(1830) relative to the K(1460). 



We have solved the Laplace transform QCD sum rules using the values: < agG"^ >~ 
0.024GeV* , < ss >^< uu >= — O.OlGeV^, and allowing Aqcd and Sq to vary in the 
range: Aqcd = 280 - 380MeV, and Sq = 4 - 8GeVl The resuhs for ms{lGeY^) are very 
stable against variations in the Laplace variable M^ over the wide range: M = 1 — 4GeV . 
In Figures 1 and 2 we show some typical results for ^^(IGeV^) as a function of M^ for, 
respectively, Aqcd = 280MeV, Sq = 4 - 6GeV^ and Aqcd = 380MeV, sq = 6 - 8GeVl 
Combining all results gives 

m,(lGeV2) = 155 ± 25MeV . (22) 

This result is consistent with the other determinations in the scalar channel, Eqs.(6)-(8). 
The error given above originates exclusively from changes in the relevant parameters, and 
does not reflect possible systematic uncertainties from the hadronic sector. 



Figure Captions 

Figure 1. The running strange quark mass rris^M^ = 1 GeV^) as a function of tlie Laplace 
variable M^, for Aqcd = 280 MeV. Upper and lower curves determine the range obtained 
by varying Sq in the interval: Sq = 4.0 — 6.0 GeV^. 

Figure 2. The running strange quark mass ms^M"^ = IGeV^) as a function of the Laplace 
variable M^, for Aqcd = 380 MeV. Upper and lower curves determine the range obtained 
by varying Sq in the interval: Sq = 6.0 — 8.0 GeV^. 
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